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Abstract —In this paper the authors describe the problem of ac¬ 
quisition of interfered signals and formulate a filtering problem. 
A frequency-selective compressed sensing technique is proposed 
as a solution to this problem. Signal acquisition is critical in 
facilitating frequency-selective compressed sensing. The authors 
propose a filtering compressed sensing parameter, which allows 
to assess if a given acquisition process makes frequency-selective 
compressed sensing possible for a given filtering problem. A 
numerical experiment which shows how the described method 
works in practice is conducted. 

Index Terms —Analog-digital conversion, Compressed sensing, 
Signal sampling, Interference reduction, 


parameter, which assesses if a used acquisition process makes 
frequency-selective compressed sensing possible. 

The paper is organized as follows. The problem considered 
in this paper is discussed in Section II. Frequency-selective 
compressed sensing is described in Section III. Filtering com¬ 
pressed sensing (fCS) parameter is proposed in the Section 
IV. The performance of the proposed solution is assessed by 
numerical simulation in Section V. The paper is concluded 
in Section VI. The paper follows the reproducible research 
paradigm [19], all the code and figures associated with the 
experiment are available online 1 . 


I. Introduction 


II. Problem formulation 


The Shannon-Nyquist sampling theorem states that perfect 
signal reconstruction of any signal requires a sampling fre¬ 
quency higher than twice the maximum frequency component 
in the signal [1]. In practical situations we thus need analog 
anti-aliasing filters prior to the analog-to-digital conversion 
(ADC) to facilitate the above [2]—[5] and to reduce the risk 
of saturating the ADC due to limited dynamic range, causing 
nonlinear distortion [4]—[6]. However, applying such analog 
filtering is design and implementation challenging [7]—[10], 
particularly in the radio frequency range. Therefore a digital 
solution would be preferred if possible. 

In this paper the authors propose a frequency-selective 
sampling method based on the compressed sensing technique. 
In recent years a new idea of signal sensing, known as com¬ 
pressed sensing (CS) has emerged [11]—[15]. This technique 
can be used to successfully reconstruct signals that are sampled 
at a sub-Nyquist rate, provided the signal is sparse in some 
domain. This technique is well elaborated, however, to the 
authors’ best knowledge there are no publications in which 
compressed sensing reconstruction selectively favors certain 
signal spectrum frequencies over others. 

Frequency-selective compressed sensing proposed in this 
paper allows for wider frequency spectrum of the sampled 
signal then just the wanted signal. Hence, relaxed input signal 
filtering is possible without increasing the sampling frequency 
to the Nyquist frequency of the input signal polluted by 
unwanted high-frequency signals. The main idea is to divide 
the compressed sensing reconstruction procedure into two 
phases: an optimization phase and a final signal reconstruction 
phase, in the latter phase a limited signal dictionary is used. 
The authors propose a filtering compressed sensing (fCS) 


Let us define a bandlimited (— B\, < f < B\,) baseband 
signal Sb £ R Mxl . The Nyquist rate of the signal Sb is /yjb = 
2I?b. The signal is polluted by an interference passband (7i h < 
/ < B x , — B x < f < —f?b) signal Si £ R Mxl . The received 
signal s x £ R Mxl is a sum of the wanted signal Sb and the 
interference signal sp s x = Sb + Si (Fig. 1). The signal s x is 
bandlimited (— B x < / < B x ), its Nyquist rate is f^ x = 2 B x . 

Due to the interference signal s h the Nyquist rate fa x of the 
received signal s x is in many applications significantly higher 
than the Nyquist rate fab of the wanted baseband signal Sb. 
To enable sampling with low frequency / s (f^tb < f s 7C fax) 
the signal s x must be filtered with a high-order low-pass filter 
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Signal poluted with unwanted interference (s x ) 

Fig. 1: Frequency spectrum of the received signal s x . The 
signal consists of the wanted signal Sb (green) and the un¬ 
wanted interference ,S| (red). The signal’s frequency range is 
(~B X <f< B x ). 






3 


(High energy dissipation, 
large chip area, design problems) 





Fig. 2: Solution 1: High-order filtering enables lowered sam¬ 
pling frequency. Large filters utilize considerable chip area 
and cause high energy dissipation and integrated circuit design 
problems. 


(High energy dissipation) 



Fig. 3: Solution 2: Low-order filtering requires high sam¬ 
pling frequency. High sampling frequency causes high energy 
dissipation and is infeasible to be implemented in certain 
applications. 
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Loosely filtered signal (si) 


Fig. 4: Frequency spectrum of the loosely filtered signal si 
which consists of the wanted signal Sb (green) and the loosely 
filtered interference signal (red) which is still partly present in 
the spectrum. A part of the unwanted frequencies was removed 
(pale red). 


which removes the unwanted interference (Fig. 2). Unfortu¬ 
nately, high-order filters cause design and integrated circuit 
implementation problems due to high energy dissipation and 
chip area required to implement these filters [8]—[10]. 

Another possibility is to ’’loosely” filter a signal with 
a low-order filter (Fig. 3). Let us consider a bandlimited 
{—B\ < / < B\) signal si, which is created by applying a 1 st - 
order filter on the received signal s x . This partly removes high- 
frequency unwanted signals, however there is still considerable 
interference content present in the filtered signal si (Fig. 4). 
The Nyquist frequency of the signal si is /(, x = 2 B\. The 


baseband B\ of the filtered signal depends on the filter’s cut¬ 
off frequency / c . The Nyquist rate /I x of the filtered signal 
si is lower than the Nyquist rate /n x of the unfiltered signal 
s x , but in many applications it is still significantly higher than 
the Nyquist rate / N b of just the wanted signal Sbi 

/m> /]I X /nx (1) 

Therefore, if a low-order filter is used, a high sampling 
frequency must be applied to the signal (Fig. 3). The high 
sampling frequency causes high energy dissipation [5]-[7] and 
may be infeasible to implement in certain applications. 

III. Frequency-Selective Compressed Sensing 

Compressed sensing is a technique which allows for signal 
sampling with frequency lower then the signal’s Nyquist rate. 
Compressed sensing is possible if the sampled signal x £ 
R Mx1 can be represented as: x = \Fv where >1/ £ R Mx2 k 
is a signal’s dictionary, v £ R 2ifxl is a sparse vector - a 
vector with only few (S) non-zero elements. The number of 
non-zero elements (S) is often called the signal’s ‘sparsity’. 
A relation between signal sparsity and compressed sensing is 
well developed in [14], 

Compressed sensing can be divided into two parts: signal 
acquisition and signal reconstruction. An observed signal y £ 
R ;Vx 1 is an outcome of the acquisition process: y = <Fx, 
where the sensing matrix £ R NxM represents the (linear) 
acquisition process. The reconstructed signal x is computed as: 
x = \&v, where v is the reconstructed sparse vector. There 
are several methods for reconstmcting the sparse vector v. 
One of the most classic is basis pursuit denoising or LASSO 
[16], [17], which is an i\ optimization process. This convex 
optimization problem can be posed as: 

v = argminllvllb+fc||y-©v ||2 (2) 

V 

The authors’ aim is to decrease the necessary sampling 
frequency in case the received signal is si with bandwidth 
B\, while only a lower-frequency part of the signal needs 
to be correctly reconstructed (Fig. 4). The filtering problem 
V(B\, -Bb, \F, S) is constituted by four parameters: baseband 
of the interfered signal B\, wanted signal baseband B^, 
dictionary v l', and signal sparsity S. 

The dictionary matrix d/ used in this paper is a Digital 
Hartley Transform (DHT) matrix [18], The matrix consists of 
2K columns indexed k £ {1,2AT}. Frequencies reflected 
by the columns of the dictionary matrix are, for the first K 
columns: /^ w = (—K + k — 1 )Sf, k £ {1,..., K}, for the last 
K columns: f£ igh = (— K + k)Sf , k £ {K + 1,..., 2AT}, where 
Sf is the frequency separation between the dictionary columns. 
The dictionary matrix used must span the full spectrum of the 
loosely sampled signal si. In a typical compressed sensing 
problem there is a need to reconstruct all the coefficients in the 
vector v correctly, but, as described in Section II, in the current 
problem it is only necessary to reconstruct the frequencies 
corresponding to the wanted signal Sb. Indices k' h £ of the 
columns of the dictionary which correspond to the signal 
Sb are within the interval: 

Kj, = {A — a + 1,..., A + a:} 


a=\B b /S f ] (3) 
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Fig. 5: The idea of frequency-selective compressed sensing 


The wanted reconstructed signal s b can be computed as: 

s b = 'S'hVh (4) 


where \P £ R Mx2K is a dictionary matrix, and <f> G R NxM j s 
a sensing matrix which represents the compressed sensing ac¬ 
quisition process. Hence, it can be stated that © = /„(•&, T'), 
where / n is a column-wise normalization function. Let us de¬ 
fine an atomic filtering compressed sensing (afCS) parameter: 


Cf = /f(©, k), 


( 6 ) 


e 


the function /f is defined in Section IV-B. The parameter Q. 
signifies how well the k th entry of the sparse vector v (4) 
will be reconstructed by a reconstruction algorithm for a given 
matrix 0 = /„(<&, \P) - the lower the parameter £®, the better 
the reconstruction of the /.: lh entry. Now let us define a filtering 
compressed sensing parameter (fCS) £® : 


£b„ = max (4") 

fceK' 


(7) 


where defined as in (3), is the set of indices of columns 
of the © matrix corresponding to the wanted signal s b . The 
filtering compressed sensing parameter £® b is calculated for a 
set of columns of 0, while the atomic filtering compressed 
sensing parameter £® is calculated for a single column of 
0. To facilitate frequency-selective compressed sensing for 
a given filtering problem V(Bl, J3 b , \fq S), one must find a 
sensing matrix $ for which the filtering compressed sensing 
(fCS) parameter ^® _ 7n(*,^) j s c ] ose to zero 


where \P b is a dictionary composed of columns k' b £ KJ, 
within the inverval as in (3). The vector v b = v[fc b ] is 
the central part of the reconstructed vector v composed by 
the elements of v corresponding to the frequency spectrum 
of the signal s b . It is clear that the reconstruction process 
(2) only needs to reconstruct the v b part of the vector v 
correctly. Therefore, the acquisition process must be tailored 
so that it brings certainty into the reconstruction of the v b 
part, while uncertainty in the reconstruction of the rest of the 
vector v is allowed (Fig. 5). Signal which is not covered by 
the dictionary is seen by a compressed sensing reconstrucion 
process as noise. Therefore, a dictionary 'T which spans all 
the spectrum of the received signal s£ must be used in the 
reconstruction process. Otherwise, the interfering signal would 
be treated by the reconstruction algorithm as noise in the 
sampled signal, which would dramatically compromise the 
quality of the reconstruction process. 

IV. Signal Acquisition for Filtering Problems 

Designing the signal acquistion part of a compressed sens¬ 
ing process is critical in facilitating frequency-selective com¬ 
pressed sensing. In this section a filtering CS parameter is 
introduced. It can be used to evaluate if a given acquisition 
process, represented by a sensing matrix <t>, is suitable for a 
filtering problem V. 

A. Evaluation of Signal Acquisition for Frequency-Selective 
Compressed Sensing 

Let us define a matrix © £ R Nx2K which consists of 
normalized columns of the matrix ©' £ R Nx2K : 

(5) 


B. Computation of The Atomic Filtering Compressed Sensing 
Parameter 

Here the authors show how to realize the function /f from 
(6) which computes an atomic filtering compressed sensing 
parameter <j® for the fc th column of the matrix 0. 

For /,: lh column Or of the matrix © let us create a projection 
matrix Pfc: 

p fc = QkifiihrHi (8) 

Let us generate a matrix fl, which contains W testing vectors 
as columns. The matrix fl £ M 2KxW is composed of normal¬ 
ized columns from the matrix fl' £ R 2KxW the elements ui' k 
of which are random Gaussian values: 

= n’(k,w), u' k w ~7V(0,1) (9) 



Sampling period [us] 

Fig. 6: Atomic filtering compressed sensing (afCS) parameters 
for columns 4 - 7 of the dictionary matrix. The parameters 
were computed for different sampling periods. 
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Let us define the matrix A £ K NxW w iji c h j s t | le product of 
the matrix © by the matrix of testing vectors Cl: 

A = 0f l (10) 

Now it is possible to compute a matrix T k £ R NxW 
contains vectors from the matrix A projected onto the k th 
column of the matrix 0 : 

r k = p k A (ii) 

The atomic filtering compressed sensing (afCS) parameter 

Cjfu, °f the fc th column, computed for the u> th testing vector is: 

Ck, w = \hkh^^(k,w)\ ( 12 ) 

where 7 ^ is the u> th column of the matrix T An estimated 
atomic filtering compressed sensing parameter of the fc th 
column of the matrix © is: 

Cf = mcLx((®J (13) 


while the the parameter computed for the second filtering 
problem V 2 is computed using atomic filtering compressed 
sensing parameters for columns {4,7}: 

CmkHz = max (Cf> Cf 1 CeP i Op) (15) 

The parameter £® Hz is close to 0 for sampling periods 25.0 //s 
and 28.5 /its (marked with red arrows in Fig. 6 ). The parameter 
CiokHz * s c l° se t° 0 for the sampling period of 25 /is (marked 
with a black arrow in Fig. 6 ). 

Average reconstruction success ratio was measured for 
filtering problems V\ and TV The ratio was measured over 
10 4 random cases. The sampling period is swept from 25 /is 
to 30 /rs with 0.5 /is step. The number A} of interfering tones 
is swept over N\ £ {0,1,...,4} for the filtering problem V\ 
and N\ £ {0,1,..., 3} for the problem TV The reconstruction 
is treated as successful if the signal-to-noise ratio of the 
reconstructed wanted signal Sb is equal to or higher than 25 
dB. The results are plotted in Fig. 7. 


It requires an infinite number (W — > 00 ) of testing wectors 
to determine the correct value of £® numerically. In practice, 
the number of testing vectors needed to determine £® with 
sufficient accuracy should be found experimentally. 

V. Numerical Experiment 

A numerical experiment was conducted to verify the idea 
practically. Loosely filtered signal sj (Fig. 2) consists of 
maximum 5 tones separated by 5 kHz, its Nyquist frequency 
is 50 kHz. The total number of tones currently present in 
the signal is not known to the recontruction algorithm. The 
signal dictionary V F used in the experiment is a discrete 
Hartley transform dictionary with 10 columns which reflect 
5 frequencies: {5 kHz, 10 kHz, ...,25 kHz}. Let us define 
two filtering problems V\ (25kHz, 5kHz, \F, N\, + N\) and 
V 2 (25kHz, 10kHz, \F. iVjj + Ni), where Ny, is the number of 
wanted tones. A} is the number of interfering tones. In the first 
problem the lowest frequency tone (5 kHz) must be correctly 
reconstructed (A}, = 1). In the second problem two frequency 
tones (5 kHz and 10 kHz) must be correctly reconstructed 
(Ab = 2). The frequency location of the A), wanted tones 
is known, while the frequency location and the number of A) 
interfering tones is not known by the reconstruction algorithm. 

The signal is uniformly sampled, to check which sampling 
period is the best for the filtering problems V\ and TV 
Atomic filtering compressed sensing (afCS) parameters ( 6 ) for 
columns 4 — 7 of the dictionary was measured for different 
sampling periods. The sampling period was swept from 20 /is 
to 30 /is with 0.5 /is step. The Nyquist frequency of the signal 
(50 kHz) corresponds to 20 /rs sampling period, so all of the 
tested uniform sampling frequencies were not higher than the 
Nyquist frequency of the loosely filtered signal s). Measured 
parameters are plotted in Fig. 6 . Columns 4 and 7 correspond 
to a 10 kHz frequency tone, columns 5 and 6 correspond to a 
5 kHz frequency tone. Filtering compressed sensing parameter 
(7) for the first filtering problem V\ is computed using 
atomic filtering compressed sensing parameters for columns 5 
and 6 : 

(14) 



The number of interfering tones 

Fig. 7: Reconstruction success ratio of filtering problem V\ 
(left) and V 2 (right) 


As expected, the reconstruction of the wanted signal from 
the filtering problem V\ for the sampling periods 25.0 /is and 
28.5 /is (red arrows) is ideal. Reconstruction of the wanted 
signal from the filtering problem V 2 is ideal for the sampling 
period 25.0 /is (black arrow). Surprisingly, no sparsity of the 
received signal is needed, since the reconstruction works well 
when the signal and interference occupy the whole spectrum. 
Instead of signal sparsity, the proposed method exploits that 
only a part of the spectrum is required to be reconstructed 
correctly. 


VI. Conclusions 

In this paper the authors have described the problem of 
acquistion of interfered signals and formulated a filtering prob¬ 
lem V. A frequency-selective compressed sensing technique 
was proposed as a solution. A filtering compressed sensing 
parameter was proposed for assessing if a given signal acqui¬ 
sition process makes frequency-selective compressed sensing 
possible for a given filtering problem. A numerical experiment 
which shows how the method works in practice was conducted. 


OfkHz — max (C?) Cf) 
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